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Abstract. We calculate radiative corrections to the flavor-changing neutral current process K — > Txt^l~ , 
both for charged and neutral kaon decays. While the soft-photon approximation is shown to work well for 
the muon channels, we discuss the necessity of further phase space cuts with electrons in the final state. 
It is also shown how to transfer our results to other decays such as ?; ^ j£^£~ or w ^ ■k'^£'^£~ . 



PACS. 13.20.Eb Decays of K mesons - 13.25.Jx Decays of other mesons 
corrections to strong- and weak-interaction processes 



13.40.Ks Electromagnetic 



1 Introduction 

The radiative non-leptonic decay K — >■ tt£~^£~ {£ = e, ^) 
offers insight into the nature of the weak interactions at 
low energies. These flavor-changing neutral current transi- 
tions are suppressed to the one-loop level in the standard 
model, where the CP-allowed channels K^ — > Tr^^+i'^ 
and Kg — >■ 'k^£'^£~ are expected to be dominated by 
one-photon exchange [K -^ t^I*)- The real-photon de- 
cay K — >■ 7r7 is forbidden by gauge invariance, and in the 
effective-theory description of chiral perturbation theory, 
also the virtual transition is suppressed to the one-loop 
level [1]. K ^ 'k£'^£~ has since been investigated in chiral 
perturbation theory beyond one- loop order [2, 3] and with 
isospin breaking included [4] . 



The least rare of these decays, K 



^±e^ 



' , the first 



to be discovered at the CERN PS [5] and subsequently 
studied with increasing precision at BNL [6], has now 
been investigated by the modern high-precision experi- 
ments E865 in Brookhaven [7] and NA48/2 at CERN [8], 
both with numbers of events of the order of 10^. For 
the muon decay channel K^ — > 7r^+/z~, a few hun- 
dred events have been recorded at BNL [9, 10] and Fer- 
milab [11], and NA48/2 is in the process of analyzing a 
sample of ^ 3000 events [12]; the Ks decays have been 
seen by NA48/1, though with a handful of events only 
so far [13, 14]. As a welcome side effect of the upcoming 
NA62 experiment dedicated to the search of the very rare 
K — ;> t:vv decays, a strong future increase in statistics for 
K — > tt£^£~ is also expected, making the availability of 
reliable (electromagnetic) radiative corrections an indis- 
pensable requirement. 

As the chiral representation of the K — ^ 777* form fac- 
tor starts at the one-loop level, a full analysis of radiative 
corrections in the framework of chiral perturbation the- 
ory seems highly impractical, as it would require a two- 
loop calculation already at leading non-trivial order. Fur- 



thermore, the counterterms appearing in such processes 
(weak non-leptonic with electromagnetic corrections at 
two loops) have never been considered, let alone the finite 
parts of such counterterms estimated phenomenologically. 
We therefore follow a more pragmatic approach akin to 
the study of soft-photon corrections in meson decays in 
Ref. [15] and restrict ourselves largely to universal radia- 
tive corrections, assuming a linear form factor; in this way, 
a technically simpler one-loop calculation suffices. Given 
that the data are still compatible with such a linear form 
factor approximation [7, 8] , this seems a reasonable proce- 
dure. Furthermore, as the only non-analyticity in the form 
factor occurring within the physical decay region is the 
7r"'"7r~ intermediate state that necessarily has to appear 
in a P-wave [2], the deviation from an analytic behavior 
is moderate, with the corresponding imaginary part rising 
slowly above the tttt threshold. 

We shall see below, after introducing the necessary for- 
malism in Sect. 2, that for the most interesting corrections 
concerning the invariant mass spectrum of the lepton- 
antilepton pair, the assumption of a specific form factor 
is actually not of high concern, as it factorizes neatly. 
These results are contained in Sect. 3. We discuss the soft- 
photon approximation for the bremsstrahlung, and how 
the appearance of coUinear singularities for the electron- 
positron final state enforces to go beyond that, introducing 
additional phase space cuts. The form factor slope only 
affects those radiative corrections that induce an ^+^~ 
asymmetry, which we investigate in Sect. 4. Despite quite 
different hadronic currents contracted with the leptonic 
ones, we demonstrate in Sect. 5 that a significant part 
of our results can be transferred directly to other meson 
decays involving lepton-antilepton pairs in the final state, 
such as the Dalitz decays 7r° — ?► 76+6^, 77 -^ ^£'^£~ , or the 
vector meson conversion decays to — >■ 7r"^+^~, (j) — > ri£~^£~ 
etc. Finally, we conclude with a summary; several of the 
more laborious formulas are relegated to the appendices. 
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2 Formalism 

We consider the decay K{k) — >• 7r(p)^+(p+)£^(p_) in 
terms of the kinematic variables s = (k—p)'^, t = (fc— p+)^, 
u = (k—p^)^. For the CP-allowed {K^ and Ks) channels, 
the K — >■ 7r7* form factor F{s) is defined as 

z / d'x e'C^-P)- {7r{p)\T{JU^)CAS=i (0)} |if (fc)) 



eF{s) 
MUM' 



[sik + p)^ ~ (k^ - p^)ik - p)^'] 



(1) 



which on the kaon and pion mass shell simplifies according 
to fc^ = M]^, p^ = M^. Here, J^!"^ is the electromagnetic 
current and Las=i the strangeness-changing non-leptonic 
weak Lagrangian. The form factors F{s) are (potentially) 
very different for K^ and Kg decays [2]; for simplicity 
we refrain from reflecting this explicitly in the notation. 
Contracting Eq. (1) with the leptonic current yields the 
decay amplitude 



M{K ^ tt£+1-) 



'Fis) 



Afi(4^) 



;u{p-){ji+ :^)v{p^ 



which, when squared and summed over fermion spins, re- 
sults in the Dalitz plot distribution 



dr = £^l^{x-.^)dsd., 



4(47r)5A// 



(3) 



K 



where a = e^/47r is the fine structure constant, A ^ 
A(M|., s, M^) with the Kahen function A(a, b, c) 



2(a6- 



ac - 



be), and v = t — u. Integrating over 



V, with the limits of the Dalitz plot given by 



< V < Vr, 



= aAi/2 , a - 



Am'^ 



(4) 

and m is the mass of the lepton involved, we find the 
spectrum with respect to the dilepton invariant mass as 



dr 

ds 



a^\F{s) 
2(47r)5M^' 



.A3/2 



all 

V 3 



(5) 



Due to the Dirac structures in the leptonic current, ra- 
diative corrections in these decays cannot easily be written 
as one overall correction factor to the Dalitz plot distri- 
bution Eq. (3) (compare Ref. [15]). This is easily seen by 
the virtual-photon exchange between lepton and antilep- 
ton, see Fig. 1, which is precisely the diagram responsible 
for the generation of the anomalous magnetic moment at 
leading order, and hence a new Dirac structure. In gen- 
eral, we write the radiative corrections to the Dalitz plot 
in the form 



dF- 



a'\F{s\ 



|2 r 



4(47r)5M^ 



(A — I' ){l+uj+uj) + i' Lj + m vuj 



dsdv , 

(6) 




K 



TT 



e+ 



^: 



K 




e+ 



Feynman graphs contributing to radiative corrections 

0,. 



Fig. 1 

in Ks — >■ Tv^i'^E through virtual- and real-photon processes 
We do not show diagrams responsible for (electromagnetic" 
wave function renormalization. 



where oj, <D, and w are even functions in v, while to is odd; 
all four of them are of order a. This leads to a corrected 
dilepton spectrum according to 



dF 
ds 



a^\Fis) 
2(47r)5M^ 



-A3/V(l-^)(l + f2), 



K 



(7) 



The relation between J7, w, and ui in Eq. (7) will be implic- 
itly assumed for all parts and approximations discussed 
for these correction factors in the following section. The 
factors a) and cu (which are only non-vanishing for the 
charged kaon decay) cancel in dF/ds upon integration of 
the terms odd in v. They do, however, induce an asym- 
metry in p (or £^i~ asymmetry) to which, in turn, w and 
u) do not contribute. This asymmetry will be defined and 
discussed in Sect. 4. Finally, we remark that, in this arti- 
cle, we do not discuss effects due to vacuum polarization 
in the one-photon exchange graphs, which can simply be 
taken care of by the properly renormalized running fine 
structure constant a. 



3 Radiative corrections for dF/ds 

3.1 K — )■ 7r/x+/2~: soft-photon approximation 

It is well known that loop corrections involving massless 
particles (e.g., photons) diverge in the infrared. These di- 
vergences are exactly canceled by the contributions from 
real-photon emission (bremsstrahlung) [16-18], which is 
typically included up to some maximum photon energy 
^max! remnants of the infrared divergence can then be 
seen by terms depending logarithmically on -Emax- The fi- 
nite, not logarithmically enhanced terms depend on the 
choice of reference frame in which this cut is applied. 
While a conventional choice in decays is the rest frame of 
the decaying particle, we would like to argue that this is 
not very advantageous in the case at hand: as the predomi- 
nant interest in this decay lies in the extraction of the form 
factor F{s) from data, a preferred option is to use a cutoff 
in the rest frame of the i^i~j system in the final state, or 
k — p = 0. In this way, the bremsstrahlung contribution — 
like the rest of the electromagnetic corrections — is essen- 
tially a function of s only and does not induce any i^£~ 
asymmetry by itself. The results we present in the follow- 
ing always refer to this choice for the photon-energy cut. 
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Emux is then given by the minimum of a cutoff energy Ecut 
and the kinematic limit, 



E, 



l|Sc 



s — Am 



max — mm-^ -t/cut, — ^ r- 

L 2y/S 



}■ 



(8) 



The bremsstrahlung calculation is typically the hard- 
est part of the radiative corrections at 0{a). It is drasti- 
cally simplified in the so-called soft-photon approxima- 
tion, which consists essentially in neglecting the pho- 
ton momentum in the overall momentum-conserving S- 
function; as a consequence, also all kinematic variables 
can be chosen identical to the non-radiative process. It is 
easy to show that the soft-photon approximation differs 
from the exact results by terms of order -Emax, i-e. repro- 
duces the logarithmic terms as well as terms of order -E^ax 
correctly. With a typical cutoff of E^ut == 20 ... 30 MeV, 
and a process such as K ^- irii'^fi^ in which all mass 
scales are significantly larger than -Emax, the soft-photon 
approximation is therefore expected to work very well. We 
shall verify this below. 

There may actually be other processes contributing to 
real-photon emission besides bremsstrahlung, i.e. the ra- 
diation of a soft photon from one of the external charged 
particles, namely so-called structure-dependent contribu- 
tions. These structure-dependent terms can be shown 
to also be of ©(Emax), hence comparable to those ne- 
glected in the soft-photon approximation. In the process 
at hand, there is a very large structure-dependent contri- 
bution from K^ — > 7r+7r° with subsequent Dalitz decay 
tt" — >■ 7e+e~, which however occurs at fixed s = M^q. To 
avoid contamination with structure-dependent photon ra- 
diation, the cutoff Emax may not be chosen too large. As 
a side remark, this should currently not be a major issue, 
as even the main (non-radiative) decay K~^ — > Tr+e^e" 
is only measured for s > M'^a to avoid background from 
K~^ — )■ tt'^tt'^ with subsequent n'^ — )■ e^e~. 

We begin by discussing the simplest of the decay chan- 
nels with respect to radiative corrections, Ks — >■ ir^ ^jl^ ii^ . 
Obviously, photons can only couple to the charged lep- 
tons in the final state. The diagrams contributing to both 
virtual- and real-photon corrections are shown in Fig. 1. 
The resulting correction factors are given by uj{Ks -^ 
/i^) = i-li|°'^_ and correspondingly for uj and i7, where 
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Fig. 2. Correction factor for dr/ds{Ks — > tv n fi~), both in 
soft-photon approximation (full line) and for the exact result 
(dashed), plus the correction factor for dr/ds{K^ — > 7r+/i+^^) 
in soft-photon approximation (dash-dotted curve). All terms 
are given for a photon-energy cutoff Ecut = 20 MeV. 



The expressions in Eq. (9) are ultraviolet finite and do not 
require a counterterm; this is a consequence of the Ward 
identity in QED. The resulting correction factor f^^^l- is 
plotted in Fig. 2. At small s, it is obviously dominated by 
the Coulomb pole in Eq. (9), 



(7^ 



al 



2(7 



(11) 



which is the leading approximation (in a) to the Gamow- 
Sommerfeld factor [19,20]. Next-to-leading order correc- 
tions (of 0(0;^)) become important as soon as a{l+a'^)/2a 
is not small any more, which, in the present case of muons 
involved, happens in the sub-keV region above threshold. 
Given the expected energy resolution of kaon decay ex- 
periments, the approximation Eq. (II) is therefore surely 
sufficient. 

The exact expression f^i+g- , without the soft-photon 
approximation, is given in Appendix A. For illustra- 
tion, we compare exact and soft-photon form for Ecut = 
20 MeV in Fig. 2: the differences are small, of the order 
of 0.2% at most. This is in accordance with a correspond- 
ing comparison for radiative corrections in iiT — >■ Stt de- 
cays [21] (which has rather similar kinematics to J^ — >■ 
7r/x~''/i^). Changing the cutoff to Ecut = 30 MeV hardly 
changes the picture. 

The generalization to the correction factor il{K^ -^ 
Tr~^ ld~^jj,~) for the dilepton spectrum is surprisingly simple. 
It is given as 



f2{K^ 



MA* 



i7, 



[2 



K-R 



(12) 



where Li(2;) is the dilogarithm or Spence function, 

log(I - t) 



U{z) 



t 



-dt 



(10) 



where fix-n is calculated exclusively from the diagrams 
displayed in Fig. 3, i.e. those graphs that would also con- 
stitute the radiative corrections in K^ -^ tt^vP, say. All 
virtual-photon corrections linking mesons and leptons, or 
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K 



yf^ 




Fig. 3. Additional Feynman graphs contributing to radiative 
corrections in dr/ds{K'^ —>■ tt'^£'^£~) through virtual (top) 
and real (bottom) photon processes. Again, diagrams con- 
tributing to wave function renormalization are not shown ex- 
plicitly. 



the interference of meson and lepton bremsstrahlung, are 
odd in i — M and hence cancel in the dilepton spectrum 
(they are part of the correction factors lu and cj to be dis- 
cussed in Sect. 4). Furthermore, f^Kir = ^k-k, as obviously 
radiative corrections only concerning the hadronic current 
do not induce a new Dirac structure, ljktt — 0. We only 
calculate luxtt in the soft-photon approximation, for which 
we find the analytic representation 



soft 
-'Kit 



^2 

TT 



1 



loa 



2p 

2E, 



log 



(1+p)^ 



^MkM^ 



Z\2 _ (1 _ p)2 

\ 



1 



^log (^ + ^)^ 
2p ^{A-pf^l 



(1+p)' 



2p °^ Z\2 - (1 - p)2 



I- s 
-2 Li 



2 Li 



2p 



2p 



P + A-1 



log 



p- A-1 

i^ + pf 
[A-pf 



-2 Li 



2 Li 



2p 



2p 



p-A + l 



p + A + 1 

-iiogr^i^ii±^^ 

1 ^A-^-il-pYMl 



(13) 



where E = (M]^ + M^)/s, A ^ (Ml - M^)/s, and 
p = X^/'^ / s. We note that also ujk-k is ultraviolet finite 
and does not require a counterterm. The contribution of 
graphs with a virtual-photon coupling to the K -^ 7r7* 
vertex, see Fig. 3, is essential for this cancelation to work. 
This vertex, with full dependence on the form factor F{s), 
can be constructed from gauge invariance requirements 
along the same lines discussed in Refs. [22, 23] (for the 
bremsstrahlung parts in Kij,^). We find therefore that all 
graphs contributing to radiative corrections in dF/ds as 
shown in Figs. 1 and 3 depend on the full form factor 
F{s) without constraints on its functional form, so F(s) 
can be fully factorized, and the correction factors il are 
independent of it. [This of course ceases to be true once 
"internal" radiative corrections are included as e.g. pho- 



ton exchange within a charged meson loop contribution to 
the form factor; see Ref. [21].] 

We also show the correction factor Q{K^ — > 7r+/i+/x^) 
(in the soft-photon approximation) in Fig. 2. The addi- 
tional term is seen to increase it compared to the neutral 
kaon channel by a smooth function contributing about 1% 
to 1.5%. 



3.2 K -^ Tve+e- 

Inserting the electron mass Trig into the correction factor 
in the soft approximation, Eq. (9), or the exact form in 
Appendix A, leads to two observations that may seem un- 
settling at first sight. First, the correction factor ile+e- is 
large, up to 10% or so, despite the fact that the Coulomb 
pole that is the dominant effect in the muon channel is 
hardly visible, as its characteristic extension is now given 
by the scale of the electron mass and hence tiny. Second, 
there are much more severe deviations between the ex- 
act form and the soft-photon approximation, even though 
both lead to results larger than naively expected. 

The explanation for this behavior can be found by ex- 
panding fle+e- from Eq. (A.l) in the electron mass, ne- 
glecting terms of 0{m), which leads to 



TT 4 



loge 



(l-e)(3-6) 



loK^ 



1-e 



(3-e)log(l-e) 

^2 



11 -3e 



-21oge+y-2Li(e)|+0(m) , (14) 

with e = 2Emax/Vs and 6 — m'^/s. The result for the 
log 5 enhanced terms is in agreement with Ref. [24]. Nu- 
merically, one finds that this reproduces the exact form to 
an accuracy of 10^^. Equation (14) contains terms that di- 
verge logarithmically in the limit of the electron mass ap- 
proaching zero, so-called Sudakov logarithms. This is not 
in contradiction with the Kinoshita-Lee-Nauenberg the- 
orem [25,26], which only states the absence of such mass 
singularities in inclusive transition probabilities. Indeed, 
it is easy to check that in the inclusive limit, correspond- 
ing to e = 1 in Eq. (14), all terms ex log (5 vanish, and the 
whole correction factor reduces to f2e+e- = a/4:TT + 0{m). 
We also understand why the soft-photon approximation 
is insufficient: it amounts to neglecting all terms of 0{e) 
in Eq. (14), which changes the coefficient of the mass log- 
arithm and fails to yield the required cancelation in the 
inclusive limit. 

The origin of the electron mass singularities is well 
understood: they stem from photons radiated collinearly 
off the (nearly massless) electrons. The close-to-on-shell 
electron propagator in a bremsstrahlung diagram, with a 
photon of momentum Z, is of the form 



1 



1 



2p^l 



m^ +p_ — |p_| cos( 



(15) 
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Fig. 4. Dalitz plot for 7* — s- e^e^7 at maximal s = (Mk — 
Mt^Y' ■ Shown are the positions of the collinear singularities for 
si = and S2 = (broad grey bands) as well as the infrared 
singularity at their intersection (thick grey dot) . Drawn relative 
to those are various cut positions in Q^±^i (dashed line), Smax 
(vertical dash-dotted), and 6^+^- (full black). 



which in the limit 



diverges in forward direction, 



cos 9 e- J = 1. To illuminate the structure of the various 
divergences, we show a sketch of the Dalitz plot for the 
decay 7*(p — fc) — )• e'^ {p+)e~ {p-)j{l) in Fig. 4, for a 
"mass" of the decaying virtual photon corresponding to 
the kinematic maximum .^Smax = Mk — M^^. We dis- 
cuss this sub-decay in terms of three Mandelstam variables 
Si/2 = (Pt +0^ S3 = {p+ +P-)^, obeying si + S2 + S3 = 
s + 2nri? « s in the massless limit, where the Dalitz plot 
degenerates into a triangle. Figure 4 shows the collinear 
singularities occurring along the lines si ss 0, S2 « 0, as 
well as the infrared singularity at Z = at the intersection 
of those two lines. It is obvious that excluding all pho- 
tons with an energy above Scut leaves mass singularities 
in the resulting corrections, stemming precisely from hard 
collinear photons. 

However, much as the inclusion of soft photons in a 
realistic measurement is an experimental requirement due 
to finite detector energy resolution, so is the inclusion 
of collinear photons due to finite angular resolution. We 
therefore have to cut not only on the photon energy, but 
also on the angles 9^±^} Figure 4 shows the position of 
such an angular cut in the e"'"e~7 Dalitz plot: obviously, 
the hard collinear photons are now part of the correction 
factor. The result for J?e+e- with such an angular cut is 



^ This of course also holds true in principle for brems- 
strahlung off heavy particles, such as the muon decay chan- 
nel discussed in the previous section. However, without any 
strong collinear enhancement, the excision of photons at rela- 
tive angles below, say, 20° from the full solid angle leads to a 
modification of the correction factor below the permille level — 
a negligible effect. 



0.00 






1 


' 1 


' 1 ' 




\^ 


0.10 


0.20 


1/2 ''•^o 
5 [GeV] ■ 




v^ 








\>-> 


V. 






-0.01 


<-.N-.\ 




- 












^^ 




v. N ^^ 






^ 




<••%. 




- 


= ■" 




v.. ^\. 






S=^ -0.02 


- 


<• 


^v ^'""--^ 


- 


t 




V. 




"^---...^.^^^ 


^ 








^^^"'^^"----^ 












t -0.03 


- 




'<- 


^^ 


a 
















'-< 




-0.04 










n n< 








^, 



Fig. 5. Correction factor for dr/ds{Ks — > 7r"e^e ) in the 
massless limit, with cuts on 9^± — 20° (full line), 10° 
(dashed), and 5° (dash-dotted), all with Scut = 20 MeV. For 
comparison, the dotted lines nearby show the same curves for 
the physical electron mass. 



given, in the massless limit (m = 0), as 



e+e 



a \ 1 



TT [4 

1-e 



(l-6)(3-6) 



loge 



(3-e)log(l-e)- 



1 — c 
log '- 



1 + c, 



■e7 



11 -3e 



1 — Cp 



1 + c, 



67 



(l-e)(l-ce^) 



log 



V 1 + Ce-y / / 



y + 2Li(e) \ , 



(16) 



with Ce^ ~ cos6g±^. Equation (16) holds for realistic 
(small) angular cuts; for cos0e±7 < ~e/(2 — e), i-e. ba- 
sically with all photons radiated in the forward direction 
relative to electron and positron cut out, the angular cut 
alone becomes more restrictive than the photon-energy 
cut, and the above formula has to be continued with the 



replacement 



'2ci=-y/{l — Cg-y). The result then trans- 



-67 



-1. 

(16) 



lates smoothly into the fully inclusive one for 
On the other hand, for small angles or c^^ « 1, Eq. 
shows how the mass singularity in Eq. (14) is translated 
into a kinematical or phase space singularity. 

/?''+ _ is displayed graphically for three different an- 
gular cuts in Fig. 5, all with i<^cut = 20 MeV. For compar- 
ison. Fig. 5 also shows, in addition to the massless result 
Eq. (16), the exact result with finite electron mass, ob- 
tained numerically: the differences are tiny. Obviously the 
correction factor is largest with the tightest cut on the an- 
gle, however for realistic angles, it is small over the whole 
dilcpton distribution range, on the level of 2% to 4%. 

There is another possibility how to apply phase space 
cuts such as to cover all the singular regions in the brems- 
strahlung that is maybe slightly less obvious, namely to 
cut on the positron-electron angle ^e+e- • All singularities. 
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lines nearby show the same curves for the physical electron 
mass. 



both the infrared divergence and the colhnear ones, occur 
for electron and positron back-to-back, i.e. cos6'e+e- ~ 
— 1. A hne of constant 9e+e- is also shown in Fig. 4. The 
result for i7e+e- with such an e^e~ colhnear cut applied, 
again in the massless limit, is given by 



^! 



TT 4 



1 



Cee + 1 
1 + Ce 



4(1 -Ce 



(1 



■log 



(17) 



with Cee = COS e+ e- ■ Equation (17) again is in accordance 
with Ref. [24]. With no cut on any dimensionful quantity 
applied (and to = 0), Q'^% _ is bound to be a constant. 
Comparison to results at finite electron mass obtained nu- 
merically for three different cut angles, see Fig. 6, show 
some deviations for very small s. Again, with realistic an- 
gular cuts, the resulting correction factors are small on 
the few-percent level. 

Finally, after this extensive discussion of fie+e- with 
various cuts and approximations, which directly yields the 
correction factor f2{Ks — >■ 7r"e''"e~), we have to remark 
on the phenomenologically most relevant case, K~^ — ^ 
7r+e+e~. As for K'^ — )■ tt+zz^/j,", we have to add the cor- 
rections on the hadronic side ^/^^ = coktv to obtain the 
full result for dF/ds. It is fully justified to use the soft- 
photon approximation Eq. (13) for this part, as it does 
not involve the lepton mass at all and hence shows no 
enhancement for ttt, — >■ 0. An analytic implementation of 
angular cuts in the e^e~7 system on f2KTr is not feasible, 
however it is also not necessary, as there are no colhnear 
enhancements, and the angular cuts only remove a very 
small part of the solid angle. Estimating the effect of a cut 
on the lepton-photon angle of 20° leads to a modification 
of the correction factor of a few times 10^* and can hence 
safely be disregarded. 
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Fig. 7. Feynman graphs contributing to the £'^£~ asymme- 
try in K+ —>■ n'^i^£~ through virtual-photon corrections. The 
real-photon contributions via interference of lepton and meson 
bremsstrahlung, see Figs. 1 and 3 respectively, are not dis- 
played again. 

4 £~^£~ asymmetry 

There is a set of radiative corrections in /I + -^ n^l^i^ 
not considered so far, namely virtual-photon corrections 
linking the hadronic and the leptonic currents, see Fig. 7, 
plus the interference terms of bremsstrahlung from leptons 
with bremsstrahlung from mesons. All these contributions 
are odd under exchange of i'^ and i~ and hence cancel in 
dr/ds, see Eq. (7). However, they induce an i~^i^ (or t—u) 
asymmetry, which we define in a normalized way as 



A.is) 



dr\ 



diysgn(h') 



dsdv 



(18) 



where sgn(a;) is the sign function. With this definition, the 
asymmetry is independent of overall coupling constants, 
and we expect it to be on the typical level of radiative 
correction effects, i.e. a few percent at most. Referring 
back to Eq. (6), we can re-express A^ in a very simple 
form in terms of the correction factors Co and w as 



A,{s) 



1 



A3/V(l - a2/3) Jo 



ay 



dv[{\- 



V ^(Jj + m vCj\ 



(19) 

The explicit results for Cj and Cj are significantly more 
complicated and lengthy than those shown so far and are 
therefore relegated to Appendix B. Several features of 
these corrections however need to be discussed in some 
detail. First, in contrast to what we found for the correc- 
tions to the dilepton spectra, w and w are not automati- 
cally ultraviolet finite. Naturally, nothing is known about 
the finite coefficients of the required counterterms, whose 
effects we therefore have to estimate by different means. 
We do this by employing the corresponding ultraviolet 
scale dependence of the loop contributions: choosing the 
mass of the p as a natural central scale, we vary it by a 
factor of e = 2.71 . . ., i.e. between 0.37Mp and 2.71Mp. 
We believe this to be a reasonable uncertainty estimate. 

Second, also in contrast to the correction factors for 
the dilepton distributions, the dependence on the form 
factor F{s) does not factorize any longer in the loop con- 
tributions of Fig. 7. In fact, a hadronic loop in F{s) would 
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immediately lead to non-trivial (i.e. non-factorizing) tvi^o- 
loop graphs for the photon corrections. As already indi- 
cated in the introduction, we approximate such effects by 
assuming a linear form factor, vi^ith the slope as measured 
in Ref. [8], F'{0)/F{0) = (2.32 ± 0.18)M^^ The cor- 
rection factors a), a) then contain terms that depend on 
F'{s)/F{s). 

There is another subtlety in the calculation of the 
graphs in Fig. 7. We have always drawn Feynman dia- 
grams with point-like Kirt^ (,~ vertices, owing to the fact 
that the first term in the form factor definition Eq. (1) 
contains a factor of s, canceling the photon propagator. 
In contrast to all graphs in Figs. 1 and 3, the second term 
in the form factor ex (fc^ — p^){k — p)^, which does not 
contain such a factor of s, does not vanish any more di- 
agram by diagram for loop contributions with photons 
linking mesons to fcrmions. The cancelation of apparent 
two-photon-cut contributions in the full result is hence a 
useful check on the calculation. 

Finally, wc only include bremsstrahlung in the soft- 
photon approximation (thereby in particular eschewing 
the non-trivial question of how to properly define the 
asymmetry, once the real-photon momentum is not ne- 
glected in the kinematic relations), also for the electron- 
positron final state. We deem this justified for two reasons: 
first, we find that the electron mass or collinear singulari- 
ties actually cancel in the asymmetry Eq. (19), hence there 
is a smooth limit of the electron mass approaching zero, 
no unnatural enhancement occurs, and there is no need 
to further cut on the angles in the phase space integral. 
Second, the remaining uncertainty inherent in the soft ap- 
proximation we consider well-covered by the sizeable (as 
we shall see below) error due to the ultraviolet scale de- 
pendence discussed above. 

The resulting asymmetries Ai,{s) arc displayed in 
Fig. 8 for K+ -> 7r+/i+^~ and in Fig. 9 for K+ -^ 



TT+e+e . Both are shown for a photon-energy cut of 
Scut = 20 McV; the variation of the cutoff by 10 MeV 
leads to changes in the asymmetries significantly below 
the ranges displayed by the grey bands, which represent 
the estimated uncertainties due to the ultraviolet coun- 
terterm. Similarly, the error due to the experimental un- 
certainty in the slope of the form factor quoted above is 
totally negligible in comparison. The main features of the 
results shown in Figs. 8 and 9 can be readily understood: 

1. The sign of the asymmetry is negative: positive v corre- 
sponds to large t and small m, or a small invariant mass 
of 7r+ and i^ , for which the repulsion of the two pos- 
itively charged particles in the final state dominates; 
vice versa, for negative v or small t, the attraction be- 
tween TT^ and £~ leads to an enhancement. 

2. The asymmetry for the muon channel vanishes for 
small s — >■ 4?7i^: due to the integration limit ^max = 
ctA^/^, the fact that lj is odd in v, plus the explicit 
factor of ly in front of w, it is obvious that the numer- 
ator in Eq. (19) has one more power of a than the 
denominator, hence A^{s) vanishes square- root-like at 
the lower end of the spectrum. This effect is invisible 
in the electron asymmetry due to the smallness of the 
electron mass, for which effectively cr « 1. 

3. A similar argument explains the behavior at the up- 
per end of the spectrum, for large s —>■ (Mk — M^r)^: 
were it only for the correction factor w, the numerator 
would be suppressed compared to the denominator by 
a factor of A^'^, and the asymmetry would approach 
zero; however the term ex w, when integrated over v, 
is enhanced by a factor of A^'^, hence the divergence 
seen in the muon asymmetry. On the other hand, uj 
comes with a prefactor of rn^, such that this term is 
negligible in the electron channel, where we therefore 
indeed see the asymmetry going to zero at the upper 
limit of phase space. 
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4. Finally, there appears to be a cusp-like structure in 
the niuon asymmetry, Fig. 8. This can be traced back 
to the V integration hitting the Coulomb poles in the 
t and u channels, which happens for a* — m(M'j^ — 
M^)/{M,r + m). The Coulomb pole is an integrable 
singularity of the form (i — tthr)^^^^ (at this order in 
a), which is touched tangentially by the integration 
(phase space) limit, leading to a non-regular behavior 
of the type 



I 



dt 



a(s) 



\/t ~ ^thr 



fc(s-S*)2 



dt 

Vt 



V(^^^ 



(20) 

hence like an absolute-value function. Again, as s* 
scales with to, this non-regular behavior practically co- 
incides with the threshold in the electron channel and 
is therefore invisible in Fig. 9. 

The overall size of the two asymmetries is indeed typi- 
cal for an electromagnetic effect on hadronic observables, 
it lies in the range of l%-2%. The asymmetry for the 
electron channel shows no enhancement, and is virtually 
identical to the one in the zero mass limit. 



5 P ^ V£+£- and V -^ P£+£- decays 

It is obvious that any decay P — )■ P'£^£^ dominated by 
one-photon exchange, where P, P' denote pseudoscalar 
mesons, will have the same universal electromagnetic cor- 
rection effects as discussed in the previous sections. They 
only depend on the structure of the hadronic current, 
Eq. (1), which in turn is dictated by parity and gauge in- 
variance, so the central results can be readily adapted for 
e.g. charm or bottom decays, with the appropriate mass 
replacements. 

The transfer of our calculations to decays of the type 
P -^ V£^£^ or y ^ P£^£^, for pseudoscalar P and vec- 
tor V mesons, may be less obvious. Physical processes 
of these types, as already mentioned in the introduction, 
include the 7r° and j] Dalitz decays tt" -^ ■ye'^e^ (see 
e.g. Ref. [27] for theoretical work on this decay), r] — > 
'y£^£^ [28], or the vector meson conversion decays such as 
w -> 7r°£+£~ or (/) -> ri£+£~ [29]. Reference [27] gives a de- 
tailed account on radiative corrections in 7r° — > je'^e~ far 
beyond the universal effects considered here, however rely- 
ing mainly on the soft approximation for the real-photon 
radiation and without an analysis of collinear singularities 
and angular phase space cuts. 

We generically denote decays of these types by A{k) -^ 
B{p)£~^{p+)£'{p-), with M^jg the masses of A and B 
{Mb = for the tt" or rj Dalitz decays). What is common 
to these channels is that they involve anomalous decay ver- 
tices of odd intrinsic parity, so the hadronic vector current 
is written in terms of the Levi-Civita tensor, 

{B{p)\J^^mA{k)) ex e^'^^PelKppFis), (21) 

where the form factor F is again a function of s = (k—p)'^, 
and we have spelled out the case of the P ^ V decay, with 



e* the polarization vector of the outgoing B. Equation (21) 
can be transformed into the V ^ P case by the simple 
replacement e* — >■ e^, e^ then being the polarization vector 
of the incoming A. 

It is conventional to normalize the rates by the real- 
photon width Pq = r{A — > B7) (which, in contrast to 
K — >■ 777, is not required to vanish by gauge invariance). 
Including radiative correction factors, we can write the 
Dalitz plot distribution as 



dP _ na\F{s)\ 
'- 



Tn 87rs(M2 



M|)3 



+ (cr^A-i/^) w dsdv , (22) 



where A = A(M^,s,M^), v — t — u defined as before, 
F{s) = F(s)/F(0) is the normalized form factor, and n = 
1 for vector meson conversion and n = 2 for pseudoscalar 
Dalitz decays, w and w are the same corrections factors of 
0{a) as discussed in various forms and approximations in 
Sect. 3. We only consider processes with uncharged bosons 
(as in all the sample processes mentioned above), so we do 
not include the correction terms odd in t — u. The resulting 
dilepton invariant mass spectrum is 



na\F{s) 



dP 



dsPo 2tts{MI- Mlf 



A-V2^(l-^)(l + r2), (23) 



where f2 is again given in terms of to and (D by the same 
relation Eq. (7) as m K ^ ■k£'^£~ . So while the Dalitz 
plot distribution Eq. (22) is entirely different from that of 
Eq. (3), owing to the different spin of one of the bosons 
involved, the dilepton mass spectrum and its radiative cor- 
rection factor look, apart from the different form factors 
involved, the same. We emphasize that this also holds for 
the bremsstrahlung contributions beyond the soft-photon 
approximation, for which this is far from obvious in the 
course of the calculation; see Appendix C for details. 



6 Summary and conclusion 

We have calculated the universal radiative corrections for 
meson decay processes with a lepton-antilepton pair in 
the final state, in particular the flavor-changing neutral 
current process K — > tt£^£^, with the following main re- 
sults. 

1. For the most interesting observable, the differential 
spectrum with respect to the invariant mass of the 
dilepton pair, radiative corrections at 0{a) are given 
as a simple sum of corrections to the hadronic and the 
leptonic current, both of which are ultraviolet finite 
without requiring a counterterm. 

2. For the muon final state, the soft-photon approxima- 
tion with a cut on the additional photon energy in the 
bremsstrahlung contribution is very accurate and leads 
to an enhancement on the l%-2% level, apart from the 
Coulomb pole close to threshold. 
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3. The (unphysical) excision of hard coUinear brems- 
strahlung photons leads to the survival of logarith- 
mic mass singularities in the electron channels. This 
requires the introduction of additional angular cuts 
in the leptonic radiative corrections, either between 
bremsstrahlung photon and electron/position, or on 
the electron-positron angle. We provide closed ana- 
lytic results for the cut dependence of these correction 
factors in the limit of the electron mass going to zero, 
which is shown to be a very good approximation. 

4. Radiative corrections in the charged decays K^ -^ 
TT~^i^£^ that link the hadronic and the leptonic cur- 
rent are odd under exchange of £'^ and £~ . They cancel 
in the dilepton spectrum, but induce an £~^£~ asym- 
metry on the level of l%-2%. This asymmetry is free 
of small mass singularities and not enhanced in the 
electron case, but requires ultraviolet renormalization 
via a counterterm that induces a sizeable uncertainty. 

5. Finally, we have shown that the universal correction 
factors for the dilepton spectra can be directly used 
also for other decays such as 7r° or 77 Dalitz decays 
or vector meson conversion decays like lu -^ Tr'^£'^£~, 
despite a different form of the hadronic vector current 
involved. 

These results ought to provide a useful and necessary in- 
gredient in future precision measurements of the decay 
processes concerned. 
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with (7e = \/{<t'^ — e)/(l — e) and e = 2E^ax/^/s- It is easy 
to see from Eq. (A.l) that Qi+e--fig+f- = 0{E^^^). The 
numerical evaluation and comparison of this result to the 
soft approximation for the muon final state is shown in 
Fig. 2. Equation (A.l) also allows to easily derive e.g. the 
fully inclusive correction factor, i.e. without phase space 
cuts on the real-photon emission, for finite mass tti; this 
corresponds to the limit e — a"^ . 



B Contributions to the £+£ asymmetry 

The radiative correction factors uj and a), defined in Eq. (6) 
and discussed numerically in Sect. 4, are given analytically 
by the expressions 
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A Correction factor f2{Ks -^ 7r°£+£-) 

In this appendix we show the analytic form of the cor- 
rection factor i^t+g- beyond the soft-photon approxima- 



tion discussed in Sect. 3.1, i.e. with the exact E-n 
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full mass dependence, but without any additional angular 
cuts. We find 
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As discussed in the main text, the K — > 777* form factor 
is approximated by a linear form, inducing the terms oc 
F'{s)/F{s) in Eq. (B.l). We have defined the following 
combinations of loop functions: 

H^^\x) ^{x-ir?- Ml) {2Ga{x) - G+(a;) + 3G:(x)] 

-2{Ml + Tr?)Gl{x) - 2{Ml - Tr?)Gl{x) - J„(x) , 
B^:f\x) = {x-m'){GU^)-G-{x)\ , 
V^ (x) ^ (m2 - m!) (m2 + m2 + x) G(i) (x) 

+ i (Ml - m! - 2s) (7m2 + 3Ml - 5x) G-{x) 



ir 
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[a.Z) Finally, the (finite) bremsstrahlung contribution given 
(in the soft-photon approximation) in terms of the one- 
Here, the tadpole loop function contains the ultraviolet parameter integral over ^ in Eq. (B.l) is specified by 
divergence in dimensional regularization. 
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[i is the running (ultraviolet) scale, and we have made 
implicit use of the notation Mg = ni for the lepton tadpole 
function An above. Of the loop functions Ja{s) and Ga{s), 
we only need the real parts in the two kinematical regimes 
< s < {Ma — to)2 or s > {Ma -I- to)2, where they are 
given by 
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As in the main text, these finite bremsstrahlung contri- 
butions are understood to be calculated for the photon- 
energy cutoff applied in the £+^"7 rest frame, or k— p — 0. 



C Details on bremsstrahlung 



In this appendix, we show some details on the calcula- 
tion of bremsstrahlung contributions to fii+t- beyond the 
with the KaUen function A = A(s, M2, m2) used here and soft-photon approximation. We concentrate on the results 
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for the correction factor to dF/ds. To describe the decay 
K{k) -^ Tr{p)i^ {p+)£~ {p-)j{l) , we make use of the ex- 
tended set of variables 
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(B.5) 
respectively, where 6{x) is the Heaviside function, and 
1 



si = {p- + If , 

{p++P-f . ti = {p + P+f 
{k-p+f , U2 = {p+p-f 

which are subject to the constraints 

s + 2n? , 
Ml 



S2 = {p+ + If , 
t2 = {k-p-f , 
(C.l) 
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Sl + S2 + S3 

s + t, +u, = Afl^ + M^+m^ + s, , i = l,2. (C.2) 
We also use the combinations of variables 
Vi=ti-Ui, i = l,2, Z' = i^i + ;^2, ^ = ^^1-2^2, (C.3) 
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functions used above are given by 
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Here, si and S2 are the chosen invariant variables for the 
integration of the Dalitz plot in the i^i^^ system (in the 
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rest frame of which the calculation is performed), and df2 
is the solid angle integral of this system relative to the vec- 
tor k = p. The numerators iV^ in Eq. (C.4), i = 1,2, 3, are 
obtained by performing the spin and polarization sums, 
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S3 



-S' 



M 



K 



Mi 



t)(. 



S2)S 



(C.5) 



Note that the non-trivial kinematic prefactors in Eq. (C.4) 
only stem from the normalization to the non-radiative de- 
cay spectrum; once this is reversed, Eqs. (C.4) and (C.5) 
contain the complete information to reconstruct the full 
kinematic dependence of the decay distribution of the 
bremsstrahlung process. 

The variables with non-trivial angular dependence are 
z^i/2 or S, A, respectively, which upon integration over 
the solid angle yield 



|^{zi, S, A\ S\ as] ^ I {Ml - M, 

(s 



S3 



{Ml - Ml 

(Sl - S2f 



32 - Si 



Si - S2 



(M' 



{Ml 



Ml 



-S3 J 
,2 A 

) +3 



(M; 



Mi 



Ml)\s 



— (ss - Am^) 
3s ^ ' '' 

A 



S3) - o('* + ^3) 



Putting pieces together, we find 
~2rr? 






where 



-/2 -2/1 -2/3 



S Ia 



h 
h 



m 
s 
1 
s 



dsids2 



(si — mi)"^ 
dsids2 



2m2 



dsids2 



(C.6) 



(C.7) 



(si 



h = ^ 



si 



m^ 



,2 



dsids2 



32 - m" 



Sl 



TTl" 



(C.8) 



Note that the integrals /i_4 are symmetric under the ex- 
change Sl ^ S2- We document the individual results for 
/i_4 for the most interesting cuts for K — ^ 7re"'"e~7, in (or 
close to) the massless limit, corresponding to the full cor- 
rections factors in Eqs. (14), (16), and (17). In all cases, 
we show the "hard" parts of the integrals, i.e. those that, 
in the combination Eq. (C.7), have to be subtracted from 
the fully inclusive correction factor fie+e- — C(/4:Tr+0{m). 



With a cut only on the photon energy, we find 
I^^"" = -I + e -log e + 0(6) , 



/h-'i = 2 Li(e) - — + 2\ogSloge + 0{S), 
/h-d = -(1 - e) [1 + log<5 - log(l - 6)] + 0{5) 



rhard 

-'4 — 



1-e 
2 

0{S) 



3e 



+ il + e){logS-\og{l-e)) 

(C.9) 



with e = 2ii'max/\/s and S = w? j s as in the main text. 
Imposing an additional cut on cos^gi^ 
the massless limit tti = 



^67) 



we find in 



/p = 0, /I'' = 2 log e log 



1 



-e7 



-2Li(e), 



-67 



iV- 


= (1 - e) 


IV- 


1-e 


2 



1 - log(l - e) - log — — ^ 

1 ~r Cg'y 



2(^-^)+l-c, 



£7 



1 +Ce7 



(l-Ce7)(l-e) 



log e 



2(1 -g) 



+ (1 + e) ( log(l - e) + lor ^ ''''' 



1 



-e7 



(C.IO) 



Finally, employing a cut on cos0g+g- ~ Cee, the hard, 
non-singular part of the bremsstrahlung integrals in the 
massless limit tti = is 



If 







Li 



1 



Cee + 1 



ree 

^4 — 



1-c, 

1 



1 



log- 

2Ce 



1 



1 



■ los 



l-h3Ce 



(C.ll) 



For completeness, we also add the equivalent of 
Eq. (C.5) for the P -> Vt^ll- and V -^ Pi+i- decays 
discussed in Sect. 5. The bremsstrahlung contribution to 
the correction factor can be written in the same form as 
Eq. (C.4) with iV, ^ Ni,i^ 1, 2, 3, A = A(M1, s, M|), 



N[ 



K 



Al m^ 



1 



Am^ 



\ + vl 



+ {s^-m'){Ml-Ml)v2\ , 

Am'^\ T 

+ A + 1^2 

s / 



aI m^ 



-{s2-m^){Ml-Ml)vA , 
N',^Ai^{s,-2m^)[l + ^)\+'l{ul + .l) 

[a + {Ml - Ml) 
(si — TV?) + (s2 — vn?) 



- 2m 1^11^2 + 7— 

2 s 



(C.12) 



12 



B. Kubis and R. Schmidt: Radiative corrections in AT — > -kH'^ i decays 



The angular integration Eq. (C.6) holds as before with the 
replacements Mk — >■ Ma^ M^^ —^ Mb- As already pointed 
out in Sect. 5, the result Eqs. (C.7) and (C.8) are then 
reproduced identically as in the decay K -^ 7r£+£~. 
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